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A novel hierarchical search technique is presented for all-sky surveys for continuous gravitational-wave sources, 
such as rapidly spinning nonaxisymmetric neutron stars. Analyzing yearlong detector data sets over realistic 
ranges of parameter space using fully coherent matched-filtering is computationally prohibitive. Thus more ef- 
ficient, so-called hierarchical techniques are essential. Traditionally, the standard hierarchical approach consists 
of dividing the data into nonoverlapping segments of which each is coherently analyzed and subsequently the 
matched-filter outputs from all segments are combined incoherently. The present work proposes to break the 
data into subsegments shorter than the desired maximum coherence time span (size of the coherence window). 
Then matched-filter outputs from the different subsegments are efficiently combined by sliding the coherence 
window in time: Subsegments whose timestamps are closer than coherence window size are combined coher- 
ently, otherwise incoherently. Compared to the standard scheme at the same coherence time baseline, data sets 
longer by about 50 — 100 % would have to be analyzed to achieve the same search sensitivity as with the sliding 
coherence window approach. Numerical simulations attest to the analytically estimated improvement. 
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I. INTRODUCTION 

Rapidly rotating neutron stars are anticipated to emit 
continuous gravitational-wave (CW) signals through various 
plausible scenarios [T'-'SI due to asymmetries. While the ma- 
jority of such neutron stars eludes electromagnetic observa- 
tions, their population might potentially be probed only by 
means of gravitational waves |7|. Currently, an international 
network of laser-interferometric detectors |8 -11J is in opera- 
tion. The observational upper limits obtained from known ra- 
dio pulsars lfT2l I13J and all-sky surveys lfT4l4T8]| akeady con- 
strain the physics of neutron stars, and thus a detection of a 
CW signal from a spinning neutron star would shed light on 
their currently rather uncertain physics 1 19 |. 

Extremely sensitive data analysis techniques are needed to 
detect prior unknown CW sources because of their expected 
low signal-to-noise ratios. A powerful method has been de- 
rived 1 20 1 based on the principle of maximum likelihood de- 
tection, leading to coherent matched filtering. CW signals 
are quasimonochromatic with a slowly changing intrinsic fre- 
quency. However, the Earth's motion relative to the solar 
system barycenter (SSB) generates a Doppler modulation in 
amplitude and phase of the waveform at a teiTestrial detec- 
tor. As shown in ll20l . the coherent matched-filtering statistic 
can be analytically maximized over the amplitude parameters 
describing the signal's amplitude variation. The so-obtained 
coherent detection statistic is referred to as the J^-statistic, 
which can also include multiple detector data streams VT\\ . 
Thus, an explicit search (evaluating F) is only done over the 
phase parameters describing the signal's phase evolution: the 
source's sky location, frequency, and frequency derivatives 
("spindowns"). 
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However, what ultimately limits the search sensitivity in 
scanning the entire sky for previously unknown CW sources 
is the finite computational resources. For yearlong data sets, 
searching a realistic portion of parameter space is computa- 
tionally absolutely impractical EOl l22l . This is due to the 
apparently enormous number of template waveforms needed 
to discretely cover the search parameter space, increasing as a 
high power of the coherent integration time. In consequence, 
viable all-sky fully coherent J^-statistic searches are restricted 
to much shorter coherent integration times, despite that the F- 
statistic can be very efficiently computed using the fast Fourier 
transform (EFT) algorithm EOl l23l . 

All-sky surveys sifting through yearlong data sets for previ- 
ously unknown isolated CW sources are accomplished by in- 
coherently combining either excess power or J^-statistic val- 
ues from shorter segments of data. In the power-combining 
methods lfl4l[T5]| , the segment duration is chosen short enough 
(typically, 30 minutes) so that the CW signal power resides in 
a single frequency bin during each segment. In contrast, so- 
called "hierarchical" J^-statistic-based methods Il24lj26l co- 
herently track the CW signal phase over longer segments (typ- 
ically of the order of a day or a few days). 

In this work, the "standard" hierarchical detection scheme 
refers to the following approach. Divide the data into nonover- 
lapping segments of duration T. Then, for a given point in 
search parameter space, the J^-statistic is computed separately 
for each segment and subsequently F values from all seg- 
ments are summed. Thus T defines the maximum time span 
of maintained phase coherence to the signal. This approach 
is efficient, because computing the coherent matched-filtering 
statistic just over T allows one to use a coarse grid of tem- 
plates in phase parameter space, compared to one required for 
the entire data set. Only when incoherently combining the 
J^-statistic results from all segments is a common ^ne grid of 
templates necessary. 

Recent progress in understanding the global correla- 



tions fZT", "28] of the J^-statistic in the phase parameters has 
lead to a substantially more sensitive hierarchical search tech- 
nique 1261 . This method has addressed a long-standing prob- 
lem, namely, the design of, and link between, the coarse and 
fine grids. A very useful geometric tool in this context is the 
concept of a metric, as first investigated in 12911301 . measur- 
ing the fractional loss in expected J^-statistic for a given sig- 
nal at a nearby grid point. While such a metric has been well 
studied for the coherent stage l22l[3n[32l . in l26l the first an- 
alytical metric for the incoherent combination step has been 
found by exploiting new coordinates on the phase parame- 
ter space. This analytical "semicoherent metric" has lately 
been further studied and extended to greater generality in 1 33 1 . 
This technology is currently also implemented and employed 
by Einstein@Home iTS], a volunteer distributed computing 
project carrying out the most sensitive all-sky CW surveys. 

The present work presents a novel hierarchical search strat- 
egy which builds on the results of |'26"33l, while further en- 
hancing the search sensitivity. Previous hierarchical search 
methods divide the data into nonoverlapping segments whose 
length is equal to coherent time baseline T. Here, a partition- 
ing of the data into segments shorter than the desired maxi- 
mum coherence length T is considered and subsequently one 
"slides" a coherence window of size T over the segments. As 
a result, segments which are closer than coherence window 
size T are coherently combined and otherwise incoherently. 
This scheme also ensures that the same semicoherent metric 
as derived in ||26][33] governs the template grid construction in 
phase parameter space but considerably enhances the overall 
search sensitivity. 

Section [Iljbriefly recaps the continuous gravitational- wave 
signal waveform. Section III describes the coherent matched- 
filtering statistic F and what we refer to as the "standard hi- 
erarchical search scheme." The idea behind the sliding coher- 



ence window approach is elucidated in Sec. IV along with an 
analytical sensitivity estimation, 
is demonstrated in Sec.JV 1 



The improved performance 
by means of Monte Carlo simula- 
compares the estimated sensitivity 
at fixed computational cost. Finally, concluding remarks fol- 
low in SeclVn] 



tions. In addition. Sec. VI 



II. CONTINUOUS GRAVITATIONAL- WAVE SIGNALS 

The dimensionless signal response function h{t) of an in- 
terferometric detector to a weak plane gravitational wave in 
the long-wavelength approximation is a linear combination of 
the form |l20l, 

h{t)=F+{t)h+{t)+Fy,{t)hy,{t). (1) 

The antenna pattern functions {t) and Fx (t) are given by 

= a(t)cos2?/; + 6(i)sin2?/;, (2a) 
Fx(i) = 6(<)cos2V' - a(i)sin2?/>, (2b) 

where ip represents the polarization angle of the signal, and 
the angle between the detector arms is assumed to be 7r/2. For 



explicit expressions of the functions a{t) and h{t), the reader 
is referred to Ref. 1201 . 

In the case of an isolated, rapidly rotating neutron star 
with a nonaxisymmetric deformation and negligible proper 
motion (cf. [34,, ,35J), the waveforms corresponding to the 
plus (+) and cross (x) polarizations are 

h+{t) = A+ sin'J'(t) , /ix W = Ax cos*(t) , (3) 

where and are the constant plus and cross polarization 
amplitude parameters, respectively, and ^'(i) is the phase of 
the signal. The parameters A+ and A-^ can be expressed in 
terms of the gravitational-wave strain tensor amplitude and 
the inclination angle i as 



A+ = /iQ (l + cos^ l) /2 , 



A, 



ho cos L . (4) 



The phase ^'(t) of the CW signal at detector time t takes the 
following form l f20l : 



= *o 



fc=0 



f{t) ■ n 
c 



k+l 



(5) 



where <i>o is the initial phase, f'^^^ = f denotes the frequency, 
and /('=>o) is the fcth frequency time derivative (also called 
"spindown"), evaluated at the SSB at reference time Iq. The 
integer s > denotes the number of frequency time deriva- 
tives to be taken into account; therefore, it holds = 0. 
The vector r{t) connects from the SSB to the detector, c is the 
speed of light, and n is a constant unit vector pointing from 
the SSB to the location of the CW source. The source's sky 
location is determined by two independent coordinates, for ex- 
ample, one can use equatorial coordinates of right ascension 
and declination, denoted by a and 5, respectively. In these 
coordinates: n = (cos 5 cos a, cos (5 sin a, sin (5). The col- 
lection of phase parameters will be summarized by the vector 
(/,/(!),. ..,/(^), a, <5). 
Using Eqs. (j2]), ([3]) and Q, it is possible to rewrite Eq. ([T]) 
as follows; 



h{t) = J2AMt) 



(6) 



where the four amplitude parameters {A-^-, Ax,ip,^o) have 
been reparametrized by the 4-vector A = (^1,^2,-^3,-44); 
whose individual components are 



Ai = A^ cos 2-0 cos<i>o — Ay. sin2-!/' sin<i>o , 
A2 — A^ sin 2'i/; cos $0 + cos 2tp sin $0 , 
A3 = —A+ cos 2%^ sin $0 ^ sin 2ip cos $0 , 
A4 = —A+ sin 2ip sin <i>o + A^ cos 2ip cos $0 , 
and the functions h^{t) have been defined as 

hi{t) = a{t) cos $(t) , h2{t) = b{t) cos <^{t) , 
hsit) = a(i)sin$(i) , hi{t) = 6(i)sin$(i) . 



(7) 



(8) 
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III. STANDARD HIERARCHICAL DETECTION SCHEME proximate ll20l . 



At a given detector time t, the detector output data time 
series is denoted by x{t). In the absence of any signal, the 
data contain only noise n{t), which is assumed to be a zero- 
mean, stationary, and Gaussian random process I.53J . When a 
signal h{t) is present, the noise is assumed to be additive, so 
that x{t) = n{t) + h{t). 

For simplicity, in this work only a single-detector input data 
stream is considered. However, based on the results of ETl 
(and also 1 32 1), it is straightforward to generalize the proposed 
search technique to multiple-detector input data, as well as to 
time- varying noise. 

To contrast with the sliding coherence window approach, 
this section describes the standard hierarchical detection 
scheme which sums one J^-statistic value from each of N 
nonoverlapping segments of duration T. For simplicity, in 
this presentation the data set is taken to be contiguous, so that 
the total data time span is written as Tdata = NT. The indi- 
vidual segments are labeled by the index j = 1, ...,N. Let tj 
denote the time midpoint of segment j, which thus spans the 
time interval [tj ~ T/2,tj + T/2]. 



A. Coherent matched-filtering of one segment 

1. The J- -statistic 

The likelihood ratio Aj for the jth segment, deciding be- 
tween the hypothesis of a signal h{t) with amplitude param- 
eters A and phase parameters p, and no signal being present, 
is written as j 



In A, = {x\h), - l{h\h)j, 



where the following inner product has been used 



(9) 



tj+T/2 



Slf Jt,-T/2 



x{t)y{t)dt. 



(10) 



with Sn^ defined as the one-sided noise spectral density for 
the jth segment. Since this work is concerned with narrow- 
bandwidth signals, S'lf ' is here taken as constant. 

As was done in Ref. EOl . the following inner products are 
combined for every segment j into a 4 x 4 matrix M whose 
components are 

M^l^ih^lK),, (11) 
where ^, = 1, 2, 3, 4. To very good accuracy, one can ap- 



ihi\h3)j « {hi\h^)j « {h2\h3)j « ih2\hi)j « 0, (12a) 



{h2\h2)j « {h4,\hi)j ~ ^Bj, 
{hi\h2)j w {h3\h4)j ~ ^Cj , 



(12b) 
(12c) 
(12d) 



with the definitions 

A, = {a\a)j, B, = {h\b),, C, = {a\h),. (13) 

In addition, we abbreviate the linear correlations {x\h^)j by 
the following compact notation: 



(14) 



Thus, Eq. (|9]) is rewritten as 



4 4 

In A, = ^ A, xbl - X E A. . (15) 
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For every segment j the log-likelihood ratio of Eq. ( 15 1 is 
analytically maximized over the amplitude parameters A. The 
maximum likelihood (ML) estimators for A obtained from the 
jth segment are denoted by 



(16) 



are explicitly given by 



J -^2 



b] 



A''' = 2 



A [j] _ 9 Bj JCg^ Cj x^f 
D., 



AlJ 

^4 



[j] 



2 _ J -^4 "-^j -^3 



(17) 



where Dj = AjBj — Cj, and Dj ^ has been assumed. 

Replacing the amplitude parameters A in In A_i of Eq. ( 15 1 
with their ML estimators A^^^ given by Eqs. ( |l7| yields the 
so-called J^-statistic for the jth segment. 



2C. 



3 ^[3] ^[3] j_ M J[3] 



^2 ^ 



(18) 



This expression can be written compactly by using the four- 
vector notation for the set of four linear correlations xji as 



J3\ = (d3] Jl3] A3] A3]\ 

A \A.-^ jAg ) 



such that Eq. ( 18 1 takes the form 



(19) 



(20) 
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where the superscript ^ indicates the transpose. Therefore, 
the J^-statistic represents a quadratic form in terms of the lin- 
ear correlations x'ji' . It should be noted that in practice Tj can 
be efficiently computed using the FFT algorithm when rewrit- 
ing the four linear correlations xji as two complex integrals; 
further details are described in Il20ll23l . 

We find that the J^-statistic can be equivalently formulated 
as a quad ratic form in terms of the ML estimators ^{f' . Using 
Eqs. \n\ to substitute the xj;?' in Eq. ( 18 1 yields 



(A1A2 + 



(21) 



which is compactly rewritten as 



is present or not, and it is exactly equal to A^I^l. It should 
also be noted that the inverse of Ai '-'l is equal to the covari- 
ance matrix of the ML estimators ^{f'. Thus, in this case, 
2J^j has noncentral distribution with 4 degrees of free- 
dom and a noncentrality parameter pj = {h\h)j, where pj 
is commonly referred to as the "optimal" signal-to-noise ratio 
(S/N). Thus, the expectation value and variance of 2J%, in this 
perfect-match case are 



Eh [2Tj] = 4 + p2 , 
where is explicitly obtained as 



4p2, (27) 



Al+Al 



' 2 
= AM^^^ A^. 



B. 



Al + Al 



Cj {A1A2 + A^Ai) 
(28) 



^ 2 



(22) 



with A representing the 4-vector of the signal's amplitude 
parameters, A — (yli,yl2j-^3,-44). Comparing Eq. (28i to 



Eq. (22 1, we find that twice the J^-statistic can be interpreted 



showing that the J'-statistic can also be viewed as quadratic ^s the ML estimator for the squared S/N: 27", 
form in terms of the ' with a coefficient matrix being equal 



to M. '■'1 . This formulation ( 22 1 of the J^-statistic is not com 



mon in the existing literature, but closely related is the work 
of 136113 81. There, the matrix is considered as a metric 
on the amplitude parameter space A, and a norm of the four- 
vector A is defined by \ \A\ \ = VXmI^M^. In this context. 



we see from Eq. ( 22 1 that the J^-statistic is simply half the 
squared "length" of the four-vector yl^^l of amplitude ML es- 
timators: Fj = ||yi[^l||V2. 



B. Incoherent combination of colierently analyzed segments 

1. The standard hierarchical detection statistic 

We denote the standard hierarchical detection statistic 
by F, which, as used in ll24l l26l [33l . represents the sum of 
one J^-statistic value Fj from each segment j. 



2. Statistical properties 



The four linear correlations xfi for a given segment j 
are Gaussian distributed random variables, whose expectation 
values and variances, respectively, are in absence of a signal, 
when x{t) = n{t), obtained as 



N N 



i=i 



(29) 



evaluated at a given fine-grid point in phase parameter space. 
Therefore, F also represents a quadratic form in terms of the 

\1\ J 
linear correlations xji , and one can compactly rewrite Eq. ( 29 1 

as 



M 



= 0, 



En 



^11 



(23) 



Therefore, in this case, the probability density function of 2Fj 
is a central distribution with 4 degrees of freedom ll20l . 
Hence, 2Fj has the following expectation value and variance, 
respectively: 



En[2Fj]^A, a. 



2J^i ,n 



(24) 



When a signal is present, which perfectly matches the tem- 
plate waveform h{t), then the expectation values correspond- 
ing to Eqs. ([23]l are obtained as 



Eh 



Eh 



- {h\h^)j , (25) 
= Mf, + m^), {h\K), . (26) 



Thus, as first noted in [20], the covariance matrix for the 
Gaussian random variables xji is the same whether a signal 



F^ ^xM 



(30) 



where the 4iV-vector x collects aU the x^^l as 



(x 



[1] J21 



(31) 



and the 47V x AN matrix Ai is defined to have the form 



M = 



V 



(32) 



Hence, the standard hierarchical scheme, using the detec- 
tion statistic F, represents the incoherent combination of 
epochwise coherent matched-filter outputs, where each epoch 
has duration T. Recall that, using Eq. ( 22 1, ^ can equivalently 



4 



be rewritten as a quadratic form in terms of the amplitude ML 
estimators A^^^ , leading to 



N 

2 ^ 



1 



■AM A' 



(33) 



where A denotes the 4A^-vector A = ( A^^^ , y^I^' , A^'^^ ). 

The problem of efficiently selecting the best coarse-grid 
value in every segment for a given fine-grid point has been 
studied in previous work 1261 [33l . and hence for rest of this pa- 
per it is assumed that such an efficient link between the coarse 
and fine grids is available. 



2. Statistical properties 

In absence of a signal, it is straightforward to show that the 
probability density function of 2^ is a central distribution 
with AN degrees of freedom ll24l[35l : hence, T has the expec- 
tation value and variance, respectively. 



2T,n 



8N. 



(34) 



On the other hand, if a signal h{t) is present, perfectly match- 
ing the template-waveform phase parameters, then ^ has a 
noncentral distribution with AN degrees of freedom and 
noncentrality parameter p^. The expectation value and vari- 
ance of respectively, are obtained as 

Eh [2J-] = m + p\ alf i^ = m + 4p2 , (35) 

where p is given by 



N 



(36) 



recalling that pj as of Eq. (28 i denotes the optimal S/N for the 
jth segment. 



It is interesting to note that p of Eq. ( 36 1 is actually equal 
to the fully coherent optimal S/N for the entire data set. The 
apparent difference in search sensitivity results from the dif- 
ferent underlying probability distributions. In the standard hi- 
erarchical scheme, there are TV times as many degrees of free- 
dom as compared to the fully coherent case. 



IV. SLIDING COHERENCE WINDOW APPROACH 

The central idea behind the sliding coherence window 
scheme is to use a window the size of the coherence time 
baseline T and to "slide" it over the data set in steps smaller 
than T to combine the coherent matched-filter outputs from 
each sliding step. This effectively amounts to the incoher- 
ent combination of coherent matched-filter outputs from over- 
lapping segments of length T, which remarkably improves 
the search sensitivity compared to the standard hierarchical 
scheme. There, only coherent matched-filter outputs from 



FIG. 1: Schematic comparison of the standard hierarchical search 
scheme and the sliding coherence window approach for the same 
coherence time baseline T. Each box represents one product of 
the linear correlations xjflip', obtained from the subsegments k 
and i, respectively. In this example, the indices take the values 
k,i = 1, 16. The light grey boxes represent the set of prod- 
ucts selected in the standard hierarchical search scheme, for A'^ = 4 
coherent segments of duration T. The dark grey boxes are the prod- 
ucts additionally selected by the sliding coherence window technique 
to enhance the search sensitivity. In this example, T is subdivided 
into g = 4 subsegments (implying a coherence overlap of 77 = 75% 
between successive sliding steps). In contrast, a fully coherent search 
over the entire data set would have to include all boxes shown. 



nonoverlapping segments of coherence length T are com- 
bined, omitting to coherently correlate large parts of the data 
which still lie within the coherence time baseline T, as illus- 
trated in Fig. [T] 

The number of templates to discretely cover the phase pa- 
rameter space searched is of great importance, since this is 
what ultimately limits the overall search sensitivity at the fi- 
nite computing power available. Therefore, it should also 
be emphasized that the sliding coherence window technique 
employs the same number of templates in phase parameter 
space as the standard scheme. The combination of overlap- 
ping coherent integrations alters neither the coherence time 
baseline T nor the total data time span Tj,,,,,. Hence, it is ob- 
vious that the same semicoherent metric as previously studied 
in 126) f33l can also be used in combination with the here- 
proposed sliding coherence window technique fS?) . 



A. Detection statistic 

For computational efficiency, the sliding coherence window 
approach subdivides every data segment of duration T into q 
"subsegments," as adumbrated in Fig.[T] Hence, each subseg- 
ment is of duration Tq = T/q, which represents the step size 
between each sliding iteration of the coherence window. 
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The subsegments are labeled by k — 1, ...,qN. Thus, in 
analogy to Eq. ( [l9] l, we define the four-vector x^'^^ for the 
fcth subsegment as xl'^l — ( jcf^' , Xj^' , Xg*^' , x|f ' ). In this con- 



text, similarly to Eq. ( [3T] i, the 4qA^-vector now collects all the 

.^^^1). Analogously, the same notation 

[fe] 



X^'^K X 

also applies to the AqN amplitude ML estimators 

As sketched in Fig.[T[ the central goal of the sliding coher- 
ence window strategy is to combine a larger number of dis- 
tinct pairs {x]'^^ , jclf' } |55 1, while still restricting the maximum 
difference between their timestamps and ti to at most the 
coherent time baseline T. Hence, the resulting selection con- 
dition is \tk - te\ < T. 

To achieve this goal, in principle an appropriate AqN x AqN 
coefficient matrix U needs to be constructed, constituting the 
following quadratic form Z, 



1 



xU x"^ , 



(37) 



which represents the detection statistic of the sliding coher- 
ence window search technique. 

In order to simplify the construction of U, we exploit the 
fact that the constants Ck are typically much smaller than the 
values of Ak and Bk, and therefore terms involving Ck are 
neglected ll56l . With this approximation, Eq. (37 i should ex- 
plicitly read as 



qN 



[tk - ti) 



1 1 3 3 



Jk]Ji] 
2 ~r -^4 -^4 



V Bk Bi 



(38) 



where the step function Qt{x) selects the pairs of linear cor- 
relations according to their time difference and the predefined 
coherent time baseline T, 




(39) 



As with J^, the detection statistic Z can also be equivalently 
reformulated as a quadratic form in terms of the amplitude ML 
estimators Aji . For practical convenience, in what follows 



we use x[f' as in Eq. (|38jl. However, note that in the above 



approximative case, xyoc ^{f ' , thus making the interchange 



between x[f' and A^^' simple if desired. 

When g = 1, it is obvious that Z coincides with the stan- 
dard hierarchical detection statistic However, if one 
chooses q > 1, the detection statistic Z is able to improve per- 
formance compared to as will be described in what follows. 

Moreover, a useful quantity is denoted by ?/, which defines 
the average "coherence overlap" between successive sliding 
steps. For the case of a contiguous data set ifSTll . as considered 
in this presentation, rj is related to q simply via rj = 1 — 1/q. 



4W 



B. Statistical properties and sensitivity estimation 

To analytically estimate the sensitivity of the sliding coher- 
ence window search, the underlying statistical properties are 



examined. Recall that, for simplicity, the data set has been 
taken as free of gaps, such that one can write the time span 
of the entire data set as Tdata = NT = qNT„. The detection 



statistic Z of Eq. (38 1 is explicitly written as 



9^ { JkY 

E 

k=l \ 

k+q-1 



AkV 



JkV 



AkV 



Ak 



Bk 



2 E 



i=k+l 



Ak]\l] ,\k]U] 



^3 ^3 



^/AkAi 



\k] m [k] u] 

2 2 4 4 

VBkBi 



(40) 



Since the actual probability density function of Z is cumber- 
some to work with, we approximate it here by a Gaussian 
distribution, which is well justified based on the generalized 
central limit theorem (provided ^ 1), as done similarly in 
previous work |25 1. Thus, we proceed by computing the mean 
and variance of Z. 

When the data consist of zero-mean stationary Gaussian 
noise only, the expectation value of 2Z is obtained as 



E,, [2Z] = 4qN , 
and the variance of 2Z is given by 



'2Z,n 



8qN{2q - 1) 



(41) 



(42) 



It is straightforward to show that for Gaussian noise, a certain 
false alarm probability PpA corresponds to a threshold Zth via 



Z,h = E„ [Z] + Gz,n V2 erfc-1 (2Pfa) , 



(43) 



where erfc denotes the complementary error function. 

Provided the presence of a signal h{t) whose phase param- 
eters perfectly match the template, then the expectation value 
of 2Z is given by 



Eh [2Z] = + p\ , 
where we defined pz as 



(44) 



Al + A" 




k+q~l 

Ak + 2 J2 ^/^kA 

e=k+i 

k+q-l 

Bk + 2 ^ V BkB, 

t=k+l 



(45) 



The probability of detection Pdet for Gaussian noise is given 
by 



2 I ^az,h 



(46) 



For current ground-based detectors, the expected CW signals 
are extremely weak, so that the small-signal situation {h <C n) 
is well justified. Thus, we approximate az.h by using az.n 
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and, by means of Eqs. (41 
the following relation: 



2 

Pz 



: £V2az,n , 
where £ has been defined as 

£ = erfc^^ (2Pfa) - erfc"^ (2Pi 



p4) i, one obtains from Eq. ( |46] l 
(47) 



(48) 



The minimum detectable gravitational-wave strain tensor am- 
plitude ho can be determined from Eq. ( 47 1, because cx 
as follows from Eq. ( 45 1. 



To obtain the estimated sensitivity scaling of the sliding co- 
herence window search in terms of the most relevant param- 
eters, the noise floor Sn is taken as constant throughout the 
data set. In addition, we replace the constants and 



by effective average values as 
and define k as 



2T, 



i A, and Bi. 



2T, 



AiAi 



Aa 



B{A2 



-Aa 



^B, 



(49) 



where the An 

are the same a,s the apart from the factor Hq 



= A^/ho. Thus, Eq. (45 i simphfies to 



p%=hlR^N{2q-l) . 



(50) 



In turn, using Eq. ( 50 1 to substitute p% in Eq. ( 47 1 and solving 
for Hq yields 



2^ 



2-' 

q 



-1/4 



(51) 



revealing the estimated sensitivity scaling of the sliding co- 



herence window search. One may further rewrite Eq. ( 5 1 1 as 
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FIG. 2: Estimated sensitivity improvement factor of the sliding 
coherence window technique over the standard hierarchical search 
strategy, shown as a function of q (number of subsegments). The 
curve is explicitly given by (2 — 1/qy^^. The horizontal dashed- 
dotted line indicates the constant value 2^^^* ~ 1.19. 



In practice, the choice of q (or equivalently rj) will generally 
have to be optimized in terms of search sensitivity at the given 
computational constraints and code implementation at hand, 
as well as for the detector data available. Further investigation 
in this direction will be presented in Sec. |VI] comparing the 
estimated search sensitivity at fixed computational cost. 



(52) 



using the previously introduced average coherence overlap rj. 



C. Comparison of sensitivity with standard scheme 



Equation (52i also reveals the estimated sensitivity im- 
provement of the sliding coherence window technique com- 
pared to the standard hierarchical search under the same as- 
sumptions. The standard hierarchical scheme is recovered for 
?7 = (i.e., q — 1). Therefore, the sliding coherence win- 
dow approach is more sensitive than the standard hierarchical 
search scheme by the factor (1 +77)^^^. In terms of q, the sen- 
sitivity improvement factor is (2 — which is shown 
inFig.|2] 

The coherence overlap r/ enhances the search sensitivity 
effectively as if increasing the number of segments in the 
standard method. In other words, to achieve the same sen- 
sitivity as with the sliding coherence window technique at 
given T, using the standard hierarchical search method, effec- 
tively rj — 1 + 1/q more segments have to be analyzed (hence 
50- 100% more data). 



V. SENSITIVITY PERFORMANCE DEMONSTRATION 

The performance improvement of the sliding coherence 
window technique is illustrated through realistic Monte Carlo 
simulations. In particular, receiver operating characteristic 
(ROC) curves are obtained to compare the standard hierarchi- 
cal search scheme that is g = 1 (corresponding to 77 = 0) and 
the sliding coherence window strategy for q = 2 (correspond- 
ing to ?7 = 50%). 

The simulated data set refers to the two LIGO 4-km detec- 
tors (HI and LI) and spans a time interval of 5 000 h. To pro- 
vide a realistic comparison, a typical value is taken for the co- 
herent time baseline of T = 50 h, which results in iV = 100. 
The software tools used are part of LALApps |f39l and em- 
ploy accurate barycentering routines with timing errors be- 
low 4/is |40|. 

In this study, the phase parameter space considered is four 
dimensional using one spindown parameter, as in current all- 
sky surveys for prior unknown CW sources ifTsHTTl . Thus a 
point in phase parameter space is labeled by p = (/,/, a, (5). 

The false alarm probabilities are found from thousands of 
different realizations of stationary Gaussian white noise with 
y/S^i = 3.25 X 10"^^ Hz"^/^. To obtain the detection prob- 
abilities, distinct CW signals with fixed gravitational-wave 
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VI. COMPARISON AT FIXED COMPUTING COST 
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FIG. 3: ROC curves comparing at fixed gravitational-wave ampli- 
tude ho the standard hierarchical search technique (lower solid curve) 
and the sliding coherence window method (upper solid curve) with 
q = 2 implying a coherence overlap of 77 = 50%. The dashed curves 
on either side of the solid curves represent estimated la errors. The 
black dotted curve is the so-called line of no discrimination. 



strain tensor amplitude of ho = 1.0 x lO^'^'' are added. The 
remaining parameters of the signal population are randomly 
drawn from uniform distributions in tp, cos l, $0, in the entire 
sky, frequencies in the interval / e [155.12, 155.16] Hz, and 
spindowns over the range of / e [—2.64, 0.264] nHz/s. 

Figure [3] compares the resulting ROC curves for the differ- 
ent search techniques. The ROC curves are computed from 
6 000 different realizations. The la errors shown in Fig.|3]are 
based on a jackknife estimate as in BTl |42I using 100 sub- 
sets. As expected, the sliding coherence window technique 
with q = 2 is substantially more "powerful" than the standard 
scheme (retrieved for q = 1), yielding a higher probability of 
detection for the same false alarm probability. 

Furthermore, the numerical results in Fig.[3]attest to the an- 



alytically estimated gain in sensitivity obtained in Eq. (51 



For example, at fixed false alarm probability of PpA = 1%, 
the achieved detection probability of the sliding coherence 
window technique is Pdet = 0.427, whereas the standard hi- 
erarchical scheme gives Pdet = 0.356. To compare these 
values with the theoretical expectation, note that Eq. ( 5 1 1 
yields VS on {2 — 1/q)^^^, where £ has been defined in 
Eq. ( [48] l and solely depends on PpA and Pdet. Thus, for 
q = 2, the predicted increase (compared to q = 1) in y/S is 
(3/2)^/'* « 10.7%. The above values obtained from the nu- 
merical simulations of Fig.[3]yield a corresponding increase in 
y/E of about 9.7%, which is in agreement with the theoretical 
prediction at the subpercent level. 



In Sec. IV B Eq. ( 5 1 1 presented the sensitivity estimate of 
the sliding coherence window technique for a given (finite) 
amount of data, disregarding aspects of computational cost 
and essentially assuming unlimited computing power avail- 
able. The present section investigates the contrary case, find- 
ing the amount of data which can be analyzed at limited 
(fixed) computational resources and a given (fixed) coherence 
time baseline T. 

The computing cost (^'^^ of a standard two-stage hier- 
archical search per a certain volume of phase parameter 
space searched can always be written as a sum in terms of 
implementation-specific constants (coh and Cincoh pertaining to 
the coherent and incoherent combination stage, respectively, 
as 



^ — 1 ^COH I5INCOH J j 



(53) 



where 7^^^ denotes the so-called refinement factor [33 1 of the 
incoherent combination stage, and A^^^^ are the number of 
segments coherently analyzed. 

The computational cost C*^"?) of the sliding coherence win- 
dow technique involves q times more summations at the inco- 
herent combination stage, thus. 



C*^) = (Co 



iNcon '7^^'^ Q 



(54) 



The total amount of data which can be analyzed at the fixed 
computational expense using the standard hierarchical search 
scheme is taken as t}:^J = T N^^\ In analogy, the amount of 
data that can be searched at given computing cost using the 
sliding coherence window approach is pflJ^ = T N'^'^\ for the 
same coherent time baseline T. 

The sensitivity of the standard hierarchical search scheme 
^^2} r)-i/4. Accordingly, the sensitivity of 



follows h'^'-' (X (T, 



scales as /i^,'^^ oc {T^' T)-V4 (2 
fine the sensitivity ratio by r 
form 



-1(9) 



the sliding coherence window technique given in Eq. (51 

l/g)-V4 



1/4 



(1) 



2- 



Thus, we de- 
which takes the 



1/4 



(55) 



At equal total computing cost, ('•'^^ — C'^^ = (, inverting 
Eqs. (53 I and (54 1 for T^} and pf^ , respectively, one obtains 



2 

q 

where the constant 



1 



1/4 



J V VI + 40 

has been defined as 

/I C CiNCOH 



'iqe-l 

1 



1/4 



2 

COH 



(56) 



(57) 



and the refinement factors have been approximated by 



(-3) 



assuming the search includes at most one spin- 
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down parameter (cf. Ref. 1331). 
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FIG. 4: Estimated sensitivity ratio r as introduced in Eq. \56\ of 
tlie standard liierarciiical searcli and the sliding colierence window 
technique at fixed total computational cost and given T, shown as a 
function of 6 defined in Eq. l[57j. The different curves correspond 
to different values of q (number of subsegments) as indicated by the 
shaded [color] bar. The dashed horizontal line corresponds to g = 1, 
for which the search methods coincide. 



Figure |4] illustrates the sensitivity ratio r of Eq. (|56]l as a 
function of and for different values of q. Two regimes are 
identified where r is slowly changing: when 9 is either very 
small or very large. This can be understood as follows. The 
specific value of depends on the code implementation, mani- 
fested in the two constants Ccoh ™d Cincoh- Thus, in the two ex- 
treme cases where one constant is much larger than the other, 
the two different limits of r result. 

First, if the implementation is such that the coherent part 
dominates the computing cost (^coh ^ Cncoh)^ this implies 
that 6 is very small. Then the sensitivity ratio is described by 



lim r — 



2- 



1/4 



1/4 



(58) 



which is the same improvement factor as given by Eq. ( 52 1. 

On the other hand, if the incoherent part is the most compu- 
tationally intensive (Cincoh ^ Ccoh)^ takes a very large value. 
In this case, the sensitivity ratio is described by 



lim r = 2 - 



1 



1/4 
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-1/8 



(59) 



It is interesting to note that, only for q — 2, the sensitivity 
ratio is always greater than 1, implying that in this case the 
sliding coherence window technique should always be more 
sensitive than the standard scheme at given T and fixed com- 
putational cost. 

In the current Einstein@Home fTS) analysis, 6 is approxi- 
mately of order 10~^. Thus, the situation in this case is rather 
comparable to the regime described by Eq. (58 1. Hence, the 



Einstein@Home search sensitivity will certainly benefit from 
employing the sliding coherence window technique. 



VII. CONCLUSION 

In summary, a novel hierarchical strategy to search for prior 
unknown continuous gravitational-wave sources has been pre- 
sented, exploiting a sliding coherence window. The standard 
hierarchical search scheme divides the data into N nonover- 
lapping segments that are coherently analyzed, and subse- 
quently matched-filter outputs are combined incoherently. 
Thereby, the duration of one segment defines the maximum 
time span of coherence. In contrast, the presented sliding co- 
herence window approach divides each of the N data seg- 
ments into q subsegments, which are thus shorter than the 
desired maximum coherence length T (size of the coherence 
window). This permits the efficient combination of matched- 
filter outputs from all subsegments in a "sliding-window" 
fashion: If subsegments are closer than T, they are combined 
coherently; otherwise, they are combined incoherently. 

As a result, the estimated search sensitivity of the sliding 
coherence window approach is considerably superior com- 
pared to the standard hierarchical scheme, while using the 
same number of coarse- and fine-giid templates to cover the 
search parameter space. At a given value of T, the sensitivity 
improvement in terms of minimum detectable gravitational- 
wave amplitude /iq scales with the fourth root of A^(2 — 1/q) 
for a contiguous data set. Since for the standard hierarchical 
method g = 1, to achieve the same sensitivity as the slid- 
ing coherence window technique between 50 — 100 % more 
data (to increase N accordingly) would have to be analyzed. 
Realistic Monte-Carlo simulations have been caiTied out con- 
firming the sensitivity enhancement. 

The sensitivity improvement can also be expressed in terms 
of the average coherence overlap 77 between successive sliding 
steps. In particular, if the data set has gaps in time, -q can be a 
useful figure of merit. The estimated sensitivity improvement 
of the sliding coherence window technique over the standard 
scheme scales as the fourth root of (1 + 77). 

In addition, the sensitivity has also been compared at fixed 
computational cost. When the computing cost of the coher- 
ent stage dominates, the above sensitivity improvement holds. 
In the case where the computational cost of the incoherent 
stage dominates, the sensitivity improvement can fade away, 
depending on the search setup. However, it is estimated that, 
when the chosen coherence window is equal to the length 
of 2 subsegments, the sensitivity is always superior to the stan- 
dard method at about the same computational expense. In gen- 
eral, the search setup (including the choice of q) will have to 
be sensitivity optimized at the given computational resources 
and the software at hand, as well as for the data available. 

However, further topics are planned to be investigated in 
future work. One of these aspects concerns the efficient im- 
plementation of the proposed technique while exploiting the 
EFT algorithm. Moreover, an optimal weighting scheme be- 
tween the coherent matched-filter outputs from different sub- 
segments could be further studied, taking into account corre- 
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lations between these. 

The sHding coherence window approach is envisioned to be 
employed by the Einstein® Home |18| project to further im- 
prove the search sensitivity of all-sky surveys for unknown 
isolated CW sources [26]. The proposed approach should 
also be extensible to CW searches for sources in binary sys- 
tems |43|. With suitable modification, the method might also 
have applicability in further related areas, for instance, re- 
garding computationally limited searches for prior unknown 
radio 1441 . X-ray Ii45il . and gamma-ray pulsars 146] - 
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